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Abstract. A brief history of the study of fluctuations in high energy nuclear collisions at the 
CERN SPS performed by NA49 is presented. The ideas and the corresponding experimental 
data on fluctuations are discussed from the point of view of their sensitivity to the onset of 
deconflnement. 



1. Introduction 

What are the phases of strongly interacting matter? This question and in particular, the 
hypothesis that at high energy densities the matter is in the form of quark-gluon plasma 
(QGP) rather than a gas of hadrons motivated the first stage of the broad experimental 
program of study of ultra-relativistic nucleus-nucleus collisions The results on the energy 
dependence of the mean hadron multiplicities and the shape of the transverse mass spectra 
obtained by the NA49 experiment at the CERN SPS indicate that the properties of the created 
matter change rapidly in the region of low SPS energies jSI (V^nn ~ 6-12 GeV). The observed 
effects confirm predictions of the transition from hadron gas to QGP ^ and thus indicate that 
in fact a new form of strongly interacting matter exists in nature at sufficiently high energy 
densities. 

Do results on energy dependence of event-by-event fluctuations confirm the conjecture of the 
onset of deconflnement at low SPS energies? This report summarizes the current status of the 
NA49 data on this subject. 

The paper starts with a brief historical review of the event-by-event physics in NA49 
(Section 2). The experimental status of the basic ideas relating the onset of deconfinement 
and the fluctuations in nuclear collisions are presented and discussed in Section 3. A summary 
given in Section 4 closes this report. 



2. A brief history of fluctuations in NA49 

The NA49 large acceptance spectrometer (Sj was constructed with the aim to study event-by- 
event fluctuations. In a letter of intend [H] it was argued: "the focus on an analysis of macroscopic 
and microscopic observables, referring to the individual events, stems from our expectation 
that the phase transition, due to the statistical fluctuations in the collision dynamics, may not 
uniformly occur in the average central Pb-|-Pb collision". The first data on central Pb-|-Pb 
collisions at 158A GeV were taken during the NA49 running period 1994-1996 and soon after 



the pioneering analysis of event-by-event fluctuations started. It appeared to be more difficult 
than originally expected. In particular the problems were caused by: 

• a need to control the experimental biases on event-by-event basis, 

• a difficulty to remove the influence of trivial fluctuations in collision geometry, 

• the absence of quantitative models predicting event-by-event fluctuations in case of the 
deconfinement phase transition. 

Numerous measures of fluctuations were proposed, see e.g. Refs. [7j, and their advantages and 
disadvantages were vividly discussed. Clearly a choice of the measure depends on the physics 
goals. It should maximize the sensitivity to the effects which we are looking for and minimize 
the sensitivity to the biasing contributions (experimental and physical). 

Due to all these problems the first results on transverse momentum fluctuations [S] and the 
fluctuations of the kaon to pion ratio [5] in central Pb+Pb collisions were published only in 
1999 and 2000, respectively. The measured fluctuations (see Figs. HandlSl) are approximately 
reproduced by a "mixed event" model in which uncorrelated production of hadrons is assumed. 
They clearly demonstrate that possible fluctuations of the early stage energy density within the 
analyzed event sample, if existent, do not lead to significant changes of the event properties. All 
collisions are "similar". 
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Figure 1. Event-by-event distribution of 
the mean transverse momentum M(pr) of 
accepted particles in the event (points) for 
central Pb-|-Pb collisions at 158A GeV. The 
solid line shows the M(pr) distribution for 
"mixed events" (histogram). 
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Figure 2. Distribution of the event-by- 
event kaon to pion ratio estimated using 
a maximum likelihood method (points) for 
central Pb+Pb collisions at 158A GeV. As a 
reference, the same procedure was applied to 
a "mixed event" sample (histogram). 



In 1997 and 1998 the data on C+C and Si+Si collisions at 158A and AOA GeV were registered. 
Together with the data on all inelastic Pb+Pb interactions at these energies, they allow to study 
system size dependence of event-by-event fluctuations. In particular, the transverse momentum 
|l()j and multiplicity fIJl fluctuations were analyzed. The data on mean pT fluctuations, see 
Fig. 131 indicate an increase of fluctuations for C+C, Si+Si and semi-peripheral Pb+Pb collisions. 
A similar effect is also observed at RHIC energies ^21- R is still lacking a well established 
interpretation. 

The hypothesis that the onset of deconflnement is located at low SPS energies motivated 
the energy scan program of the NA49 CoUaboration ^Ei- Within this program data on central 
Pb+Pb collisions at 20A, 30A, AOA and 80A GeV were registered in runs in 1999, 2000 and 2002. 
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Figure 3. The system size dependence of transverse momentum fluctuations at 158^4 GeV. Tlie 
points correspond to the NA49 experimental data 'OT, whereas the solid line to the predictions 
of the string-hadronic model Hijing j26j . 

The predicted (4j rapid changes in the energy dependence of hadron properties averaged over 
central Pb+Pb collisions, like pion and kaon mean multiplicities and transverse mass spectra, 
were observed (Hj- The search for the expected anomalous behavior of other quantities is in 
progress. Among them event-by-event fluctuations are of particular importance. This is because 
they are expected to be sensitive to the physics properties of the system which are not accessible 
in the study of the mean quantities. In the following section we review current status of the 
ideas and the experimental results on the energy dependence of electric charge, multiplicity and 
strangeness fluctuations in the domain of the onset of deconfinement. 

3. Fluctuations and the onset of deconfinement 

3.1. Net-electric charge fluctuations 

A simple estimate of the net electric charge fluctuations show that they are much smaller in 
a Quark-Gluon Plasma than in a hadron gas This difference is caused by smaller charge 
units in the QGP (fractional charges of quarks) than in the hadron gas. Thus a decrease of 
the net-charge fluctuations is expected when the collision energy crosses the threshold for the 
deconfinement phase transition. A schematic sketch of this naive prediction is presented in 
Fig.ia 

The experimental results on central Pb-f-Pb collisions at 20^, 30^4, 40A, 80^ and 158^ 
GeV shown in Fig. [3 suggest only a very weak, if any, dependence on energy. They are close 
to those expected for a gas of pions correlated only by global charge conservation (A<^q = 0) 
and significantly above the naive prediction for the QGP A<I>q ~ —0.5 |17j . These results 
do not necessarily exclude reduced fluctuations in the QGP because these can be masked by 
contributions from resonance decays. In fact a simple model of the QGP hadronization and 
resonance decay can quantitatively explain the magnitude of the measured fluctuations, see 
Fig. El 
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Figure 4. A naive expectation for the energy dependence of the net-electric charge fluctuations 
in the region of the onset of deconfinement. 
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Figure 5. The energy dependence of the net-charge fluctuations measured by A$q |17j in 
central Pb-|-Pb collisions for two different rapidity intervals Ay = 1.2 (left) and Ay = 3 (right). 
For more details see the original publication |18j . 



The influence of resonance decays on charge fluctuations depends on the size of the rapidity 
interval, Ay, in which fluctuations are calculated. If Ay is much bigger than the typical distance 
in rapidity of the daughter particles, AygMEARi the charge within the interval will not be changed 
by the decay and therefore the charge fluctuations should not be affected. On the other hand, if 
Ay is small, a large fraction of daughter particles will leave the interval and the initial net-charge 
will be significantly changed. The mean rapidity difference of two pions originating from decays 
of p{770) meson is approximately one unit of rapidity. Therefore in order to minimize the decay 
effect Ay should be much larger than 1. However, this constraint is difficult to fulfill at SPS and 
lower energies because the width of the rapidity distribution of all produced particles, AyxoT is 
not much broader than 1. Hence a rapidity interval which is large enough to be unaffected by the 
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Figure 6. The dependence of A^q on the fraction of accepted particles in central Pb+Pb 
collisions at 20-158A GeV. The prediction for the ideal QGP is indicated by the dashed curve 
(QGP), whereas the prediction for the QGP including hadronization and resonance decay is 
shown by the dotted curve (QGP+hadronization). For more details see the original publication 
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resonance decays would contain almost all particles produced in a collision. The net-charge in 
this interval would then reflect the number of participant protons and the fluctuations would be 
determined by fluctuations of the collision centrality and not the particle production mechanism. 
A sketch which illustrates these considerations is shown in Fig. [Tj At the SPS and AGS energies 
Ay w AysM ear ~ AyxoT ~ 1 and consequently the net-charge fluctuations are not sensitive to 
the initial QGP fluctuations. Hence their energy dependence may be not affected by the onset 
of deconfinement. 

3.2. Multiplicity fluctuations 

It was suggested that properly filtered multiplicity fluctuations should be sensitive to the 
equation of state at the early stage of the collision and thus to its changes in the deconfinement 
phase transition region This idea follows from the expectation that the early stage energy 
density changes from collision to collision and thus it causes fluctuations of the thermodynamical 
parameters of the matter like temperature, pressure, entropy and strangeness. Obviously the 
relation between energy density fluctuations and fluctuations of the other thermodynamical 
parameters depends on the equation of state. This opens a possibility to study equation of state 
via properly defined analysis of event-by-event fluctuations. This basic idea is sketched in Fig. |H1 
for the case of the temperature fluctuations in the phase transition domain. 

As the entropy is closely related to the particle multiplicity, and it is expected to be 
approximately conserved during the evolution of the matter created at the early stage, the 
fluctuations of entropy are of primary interest. Based on the second principle of thermodynamics 
it can be shown that the relative entropy fluctuations are given by: 
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Figure 7. A sketch illustrating different rapidity scales important in the study of the net-charge 
fluctuations as a possible signal of QGP. Ay is the acceptance window in rapidity, AysMEAR is 
a typical distance in rapidity of decay products of resonances and AyxoT is the width of the 
rapidity distribution of all particles. The observation of the predicted suppression of the net- 
charge fluctuations in QGP requires: AysMEAR < Ay < AyxoT- This condition is not fulfilled 
at the SPS energies where AysMEAR ~ Ay « AyxoT ~ 1- 
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Figure 8. A sketch illustrating how the early stage energy density fluctuations lead to 
temperature fluctuations which depend on the equation of state. The solid line shows equation 
of state with the 1st order phase transition. The width of the vertical bands indicates early 
stage energy density fluctuations for three different collision energies, whereas the width of the 
horizontal bands gives the corresponding temperature fluctuations. 



where S and E denote entropy and energy, respectively, p is pressure and e energy density. For 
the experimental study the measure of multiplicity fluctuations, Re was proposed jl9j : 



where 5n and 6E are dynamical fluctuations of particle multiplicity and energy, respectively. 
They are related to the variance of particle multiplicity V{n) as: 

y(n) = {6nf + {i6nf), (3) 

where (((5n)^) are dynamically averaged statistical fluctuations which can be measured using the 
sub-event method [SHI- The measure Re is directly related to the early stage equation of state 

i^e = (l + ^)-^ (4) 

and therefore it should be sensitive to the onset of deconfinement. As it is defined as a ratio of 
relative multiplicity and energy fluctuations it is insensitive to the magnitude of the early stage 
energy density fluctuations. One should underline that the method assumes that the early stage 
entropy fluctuations are proportional to the dynamical multiplicity fluctuations, but it does not 
require knowledge of the value of the proportionality factor. This is a clear advantage in regard 
to the methods of studing the equation of state based on the analysis of the average quantities. 

The energy dependence of Re calculated with the Statistical Model of the Early Stage 
(SMES) [4, is presented in Fig. |^ It shows a characteristic 'shark fin' structure caused 




Figure 9. The energy dependence {F = Figure 10. The scaled variance for 
{^/s — 2m)'^^^ / s^/^) of the properly filtered positively (+) and negatively (-) charged 
multiphcity fluctuations [Hj, Re, calculated particles calculated within canonical and 
with SMES [HUH]. The 'shark fln' structure grand canonical ensembles as a function of 
at the low SPS energies is caused by the onset the mean multiplicity The results 

of deconflnement. corresponds to a system with net charge equal 

to zero. 

by large fluctuations in the mixed phase region. 

This simple idea motivated further theoretical and experimental effort directed toward 
understanding of multiplicity fluctuations in nucleus-nucleus collisions at high energies. 

The effect of conservation laws on the multiplicity fluctuations were studied within canonical 
and micro-canonical ensembles. This led to an unexpected observation that the scaled 



variance of the multiplicity distribution, w = V{n)/{n), is different in various statistical 
ensembles (grand canonical, canonical and micro-canonical) even in the thermodynamical limit, 
y — > oo. The simplest example is presented in Fig. Ilfl( where the scaled variance for positively 
and negatively charged particles calculated within grand canonical and canonical ensembles 
is plotted as a function of mean particle multiplicity j2j. A detailed study of multiplicities 
fluctuations in different statistical ensembles is in progress |22| I24j . 

An important contribution to multiplicity fluctuations may come from the fluctuations in the 
collision geometry. Clearly an interaction with a small number of participating nucleons will 
result in much lower particle multiplicity than a central collision in which almost all nucleons are 
involved in the reaction process. Thus in order to extract the dynamical effects a precise control 
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Figure 11. The scaled variance of the multiplicity distribution for the negatively charged 
hadrons produced in Pb+Pb collisions (circles) and p+p interactions (triangle) at 158A GeV 
as a function of the number of participants nucleons from the projectile The solid line 

indicates predictions of a string-hadronic model Hijing j26j . 

of the geometrical fluctuations is needed. In NA49 this is achieved by selecting events with 
fixed number of projectile participants, Np^^"^ , measured by the forward VETO calorimeter 
|25j . The presented results are corrected for the resolution of the VETO calorimeter and the 
finite bin size in the measured forward energy. In contrary to the naive expectations |2S] it was 
found that the scaled variance of the multiplicity distribution increases with a decreasing 
centrality of the Pb+Pb collisions at 158^ GeV. The preliminary results^ are shown in Fig. 1111 
The observed non-trivial centrality dependence focused the efforts on its tests and understanding 
and consequently the analysis of the energy dependence of multiplicity fluctuations was delayed. 

3. 3. Strangeness fluctuations 

Similar to entropy (multiplicity) fluctuations, strangeness fluctuations are also sensitive to the 
early stage energy density fluctuations and via these to the equation of state at the early stage 
of the collisions [^. The measure -R^/g proposed to study them is defined as: 

^ Note that an error was found in the analysis of muhipHcity fluctuations in C+C and Si+Si collisions at 
158^ GeV, and thus the corresponding preliminary results shown in Refs. I14| should not be used any 

more. 
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Figure 12. The energy dependence [F = 
(■^s — 2m)^/^/s"^/^) of the properly filtered 
strangeness fluctuations, -Rg/e, calculated 
with SMES HHIj. The 'tooth' structure at 
the low SPS energies is caused by by a change 
of fluctuations in the mixed phase region. 



Figure 13. The energy dependence of the 
dynamical strangeness fluctuations in central 
Pb+Pb collisions at the SPS energies [2H] and 
central Au+Au collisions at the top RHIC 
energy [21 



where Suk and 5n are dynamical fluctuations of the multiplicity of kaons and pions, respectively. 
The energy dependence of the measure Rg/e calculated within the SMES model [H EZl is shown 
in Fig. 1121 The observed rapid change of the behavior is caused by the onset of deconfinement. 
The first results of NA49 [2Sj and STAR [22] concerning strangeness fluctuations in central 
Pb+Pb collisions at the SPS energies are presented in Fig. El The data seem to also suggest a 
change of energy dependence in SPS energies. However, a direct comparison between the data 
and the model predictions is not possible as the used experimental measure of fluctuations is 
different than Rg/e- 



4. Summary 

The basic ideas concerning fluctuation signals of the onset of deconflnement at the SPS energies 
were reviewed. The status of their experimental tests based on the NA49 data was presented. 

The naively expected signiflcant suppression of the net-charge fluctuations is not observed. 
This may be due to dilution of the initial charge fluctuations by resonance decays. 

The study of energy dependence of multiplicity fluctuations is in progress. It focuses on the 
search for the "shark fin" structure expected at the onset of deconfinement. In this respect 
two unexpected observations were made recently. Theory: the scaled variance of multiplicity 
distribution calculated within various statistical ensembles (micro-canonical, canonical and 
grand canonical) was found to be different even in the thermodynamical limit. Experiment: 
the scaled variance of multiplicity distribution for Pb+Pb collisions at 158^4 GeV increases 
significantly when going from central to semi-peripheral collisions. 



The preliminary results on the dynamical fluctuations of the kaon to pion ratio in central 
Pb+Pb collisions indicate a change of behavior at the SPS energies. The relation of this 
interesting effect to the onset of deconfinement is still unclear. This is because a proper 
comparison between the data and the model predictions is still missing. 
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